Observation of Direct Photons in Central 158 A GeV 208p b+ 2os pb Collisions 
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A measurement of direct photon production in 208 Pb+ 208 Pb collisions at 158 A GeV has been carried out in 
the CERN WA98 experiment. The invariant yield of direct photons in central collisions is extracted as a function 
of transverse momentum in the interval 0.5 < pr < 4 GeV/c. A significant direct photon signal, compared 
to statistical and systematical errors, is seen at pr > 1.5 GeV/c. The results constitute the first observation of 
direct photons in ultrarelativistic heavy-ion collisions which could be significant for diagnosis of quark gluon 
plasma formation. 
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The observation of a new phase of strongly interacting mat- 
ter, the quark gluon plasma (QGP), is one of the most im- 
portant goals of current nuclear physics research. An exten- 
sive experimental program has been undertaken at the CERN 
SPS accelerator with Pb-ion beams of 158AGeV to search 
for and investigate the QGP. Several observations, such as 
suppression of the J/ij) resonance ^ and the enhancement 
of strangeness [0], hint at an interesting new behavior of the 
matter produced in these collisions. While such observations 
imply a hot and dense initial phase with strong rescattering, 
consistent with the assumption that a quark gluon plasma was 
formed, a direct signature of the plasma and its properties is 
missing. It is therefore of great interest to search for photons 
emitted directly from the early hot phase of the relativistic 
heavy-ion collisions. 

Photons (both real and virtual) were one of the earliest pro- 
posed signatures for QGP formation [^|^] . Real photons are 
dominantly produced by scatterings of charged particles dur- 
ing the collision. Once produced, they interact with the sur- 
rounding matter only electromagnetically resulting in a long 
mean free path. They are therefore likely to escape from 
the system directly after production without further interac- 
tion, unlike hadrons. Thus, photons carry information on their 
emitting sources from throughout the entire collision history, 
including the initial hot and dense phase. 

Following early estimates of photon emission rates 
Kapusta et al. |J made detailed comparisons of the emissivity 
of the QGP and a hadron gas as two contrasting scenarios. It 
was demonstrated that the thermal emission rates of a hadron 
gas and a QGP were very similar and dependent essentially 
only on the temperature T. This led the authors to conclude 
that direct photons are a good thermometer for strongly inter- 
acting matter, but would not in themselves allow to distinguish 
between the two scenarios. 

Recently, it was shown by Aurenche et al. [|Io|] that pho- 
ton production rates in the QGP when calculated up to two 
loop diagrams, are considerably greater than the earlier low- 
est order estimates. A new higher order process of qq anni- 
hilation with rescattering was found to dominate the photon 
emission rate from quark matter at high photon energies. Fol- 
lowing this result, Srivastava [ |TTj ] has reinvestigated the pre- 
dicted photon production in heavy-ion collisions and shown 
that at sufficiently high initial temperatures the photon yield 
from quark matter may significantly exceed the contribution 
from the hadronic matter to provide a direct probe of the quark 
matter phase. 

A large number of measurements of prompt photon pro- 
duction at high transverse momentum (pr > 3 GeV/c) exist 
for proton-proton, proton-antiproton, and proton-nucleus col- 
lisions (see e.g. jl2t|). To a great extent, especially at higher 
y^i, these data can be successfully described by perturbative 
QCD calculations and provide an important foundation from 
which to study photon production in nucleus-nucleus colli- 
sions. First attempts to observe direct photon production in 
ultrarelativistic heavy-ion collisions with oxygen and sulphur 
beams found no significant excess [jO-ttq]. The WA80 collab- 



oration [ |l6| ] provided the most interesting result with a px de- 
pendent upper limit on the direct photon production in S+Au 
collisions at 200AGeV. This result was subsequently used by 
several authors to rule out a simple version of the hadron gas 
scenario [|l7|-|2"o|| and has been interpreted to set an upper limit 
on the initial temperature of T % = 250 MeV [ElJ. 

In this paper we report on the first observation of direct pho- 
ton production in ultrarelativistic heavy-ion collisions. The 
results are from the CERN experiment WA98 p^ ] which con- 
sists of large acceptance photon and hadron spectrometers. 
In addition, several other large acceptance devices allow to 
measure various global variables on an event-by-event basis 
for event characterization. Photons are measured with the 
WA98 lead-glass photon detector, LEDA, which consisted of 
10,080 individual modules with photomultiplier readout. The 
detector was located at a distance of 21.5 m from the target 
and covered the pseudorapidity interval 2.35 < rj < 2.95 
{Hem — 2.9). The particle identification was supplemented by 
a charged particle veto detector in front of LEDA. 

The results presented here were obtained from an analysis 
of the data taken with Pb beams in 1995 and 1996. The 20% 
most peripheral and the 10% most central reactions have been 
selected from the minimum bias cross section (a m in.Mas ~ 
6300 mb) using the measured transverse energy Et- In total, 
« 6.7 • 10 6 central and w 4.3 • 10 6 peripheral reactions have 
been analyzed. 

The extraction of direct photons in the high multiplicity en- 
vironment of heavy-ion collisions must be performed on a sta- 
tistical basis by comparison of the measured inclusive photon 
spectra to the background expected from hadronic decays. In- 
dividual photons cannot be tagged as isolated direct photons 
in these reactions due to the high multiplicities. To obtain 
the direct photon spectrum the following steps are performed 
(for a detailed description of the detectors and the analysis 
procedure see [^3|]): First, the raw photon spectra are accu- 
mulated after application of the photon identification criteria 
(such as transverse shower size) to the showers observed in the 
LEDA. The raw photon spectra are then corrected for contam- 
ination by charged and neutral hadrons, for conversions, for 
the identification efficiency, and acceptance. The efficiency 
includes all effects of the detector response such as distor- 
tions by shower overlap, dead and bad modules, and energy 
resolution. Neutral pions are reconstructed via their 77 decay 
branch. Invariant mass spectra are accumulated for all pho- 
ton pairs for each pair px bin. The photon-pair combinatorial 
background is estimated by event-mixing and then subtracted 
from the real-pair spectra. The yield in the 7r° mass peak is 
extracted to obtain the raw neutral pion pr spectra. These are 
then corrected for conversions, for the 7r° identification effi- 
ciency, and for geometrical acceptance. In addition, r\ mesons 
are extracted in a limited transverse momentum range with an 
analogous procedure. 

The final measured inclusive photon spectra are then com- 
pared to the calculated background photon spectra to check 
for a possible photon excess beyond that from long-lived ra- 
diative decays. The background calculation is based on the 
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FIG. 1. The inclusive photon (circles) and n° (squares) trans- 
verse momentum distributions for peripheral (open points) and cen- 
tral (solid points) 158 A GeV 208 Pb+ 208 Pb collisions. The data have 
been corrected for efficiency and acceptance. Only statistical errors 
are shown. 



FIG. 2. The 7Meas/7Bk g d ratio as a function of transverse mo- 
mentum for peripheral (part a)) and central (part b)) 158 A GeV 
208 Pb+ 208 Pb collisions. The errors on the data points indicate the 
statistical errors only. The pr -dependent systematical errors are in- 
dicated by the shaded bands. 



measured ir° spectra and the measured 77/7r°-ratio. The spec- 
tral shapes of other hadrons having radiative decays are cal- 
culated assuming my-scaling [|^| with yields relative to 7T°'s 
taken from the literature. It should be noted that the measured 
contribution (from 7r° and rf) amounts to w 97% of the total 
photon background. 

Fig. [j] shows the fully corrected inclusive photon spectra 
for peripheral and central collisions. The spectra cover the pt 
range of 0.3 — 4.0GeV/c (slightly less for peripheral colli- 
sions) and extend over six orders of magnitude. Fig. [j] also 
shows the distributions of neutral pions which extend over a 
similar momentum range with slightly larger statistical errors. 

The ratio of measured photons to calculated background 
photons is displayed in Fig. || as a function of transverse mo- 
mentum. The upper plot shows the ratio for peripheral colli- 
sions which is seen to be compatible with one, i.e. no indi- 
cation of a direct photon excess is observed. The lower plot 
shows the same ratio for central collisions. It rises from a 
value of w 1 at low px to exhibit an excess of about 20% at 
high p T . 

A careful study of possible systematical errors is crucial 
for the direct photon analysis. The various sources of sys- 
tematical errors have been investigated and are summarized 
in Table |. The largest contributions are from the 7 and ir° 



identification efficiencies and the uncertainties related to the 
77 measurement. It should be emphasized that the inclusive 
photon and neutral meson (the basis for the background calcu- 
lation) yields have been extracted from the same detector for 
exactly the same data sample. This decreases the sensitivity 
to many detector related errors and eliminates all errors asso- 
ciated with trigger bias or absolute yield normalization. The 
estimate of the systematical errors has been checked by per- 
forming the entire analysis with various photon selection cri- 
teria which change the efficiency and background corrections 
by factors of 2-3. The final results were verified to be con- 
sistent within the systematical errors for the different analysis 
cuts. Full details on the systematical error estimates are given 
in [|23|]. The total py-dependent systematical errors are shown 
by the shaded regions in Fig. |^. A significant photon excess is 
clearly observed in central collisions forp-r > 1.5 GeV/c. 

The final invariant direct photon yield per central collision 
is presented in Fig. Bl The statistical and asymmetric system- 
atical errors of FigTy are added in quadrature to obtain the 
total upper and lower errors shown in Fig. || An additional 
Pt -dependent error is included to account for that portion of 
the uncertainty in the energy scale which cancels in the ratios. 
In the case that the lower error is less than zero a downward 
arrow is shown with the tail of the arrow indicating the 90% 



3 



confidence level upper limit (j Excess + 1-28 <7jj vper ). 

No published prompt photon results exist for proton- 
induced reactions at the y/s of the present measurement. In- 
stead, prompt photon yields for proton-induced reactions on 
fixed targets at 200 GeV are shown in Fig. |] for compari- 
son. Results are shown from FNAL experiment E704 [ p5| ] 
for proton-proton reactions, and from FNAL experiment 
E629 [||] and CERN SPS experiment NA3 for proton- 
carbon reactions. These results have been divided by the to- 
tal pp inelastic cross section (<Ji n t — 30 mb) and by the mass 
number of the target to obtain the invariant direct photon yield 
per nucleon-nucleon collision. They have then been multi- 
plied by the calculated average number of nucleon-nucleon 
collisions (660) for the central Pb+Pb event selection for com- 
parison with the present measurements. This scaling is esti- 
mated to have an uncertainty of less than 10%. The proton- 
induced results have also been scaled from yfs — 19.4 GeV 
to the lower y/s — 17.3 GeV of the present measurement 
under the assumption that Ed^a^/dp 3 — J(xt)/s 2 , where 



xt = 2"pT I \fs [B8h. The y^-scaling effectively reduces the 
19.4 GeV proton-induced results by about a factor of two. 
This comparison indicates that the observed direct photon pro- 
duction in central 208 pb+ 208 Pb collisions has a shape similar 
to that expected for proton-induced reactions at the same y/s 
but a yield which is enhanced. 

In summary, the first observation of direct photons in ultra- 
relativistic heavy-ion collisions has been presented. While pe- 
ripheral Pb+Pb collisions exhibit no significant photon excess, 
the 10% most central reactions show a clear excess of direct 
photons in the range of p? greater than about 1.5 GeV/c. The 
invariant direct photon multiplicity as a function of transverse 
momentum was presented for central 208 pb+ 208 Pb collisions 
and compared to proton-induced results at similar incident en- 
ergy. The comparison suggests excess direct photon produc- 
tion in central 208 pb+ 208 Pb collisions beyond that expected 
from proton-induced reactions. The result suggests modifi- 
cation of the prompt photon production in nucleus-nucleus 
collisions, or additional contributions from pre-equilibrium or 
thermal photon emission. The result should provide a strin- 
gent test for different reaction scenarios, including those with 
quark gluon plasma formation, and may provide information 
on the initial temperature attained in these collisions. 
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FIG. 3. The invariant direct photon multiplicity for central 
158 A GeV 208 Pb+ 208 Pb collisions. The error bars indicate the com- 
bined statistical and systematical errors. Data points with downward 
arrows indicate unbounded 90% CL upper limits. Results of several 
direct photon measurements for proton-induced reactions have been 
scaled to central 208 Pb+ 208 Pb collisions for comparison. 



the INTAS under Contract INTAS-97-0158, ORISE, Grant- 
in-Aid for Scientific Research (Specially Promoted Research 
& International Scientific Research) of the Ministry of Edu- 
cation, Science and Culture, the University of Tsukuba Spe- 
cial Research Projects, and the JSPS Research Fellowships for 
Young Scientists. ORNL is managed by UT-Battelle, LLC, 
for the U.S. Department of Energy under contract DE-AC05- 
00OR22725. The MIT group has been supported by the US 
Dept. of Energy under the cooperative agreement DE-FC02- 
94ER40818. 



[1] M.C. Abreu et al., Phys. Lett. B 410, 337 (1997). 

[2] E. Andersen et al., Phys. Lett. B 449, 401 (1999). 

[3] EL. Feinberg, Nuovo Cimento 34 A, 391 (1976). 

[4] E. Shuryak, Phys. Lett. B 78, 150 (1978). 

[5] K. Kajantie and H.I. Miettinen, Z. Phys. C 9, 341 (1981). 

[6] F. Halzen and H.C. Liu, Phys. Rev. D 25, 1842 (1982). 

[7] B. Sinha, Phys. Lett. B 128, 91 (1983). 

[8] L.D. McLerran and T. Toimela, Phys. Rev. D 31, 545 (1985). 

[9] J. Kapusta, P. Lichard, and D. Seibert, Phys. Rev. D 44, 2774 



4 



TABLE I. Various sources of systematical error in the WA98 158 A GeV Pb+ Pb direct photon analysis specified as a percentage of 
(7/7r°)Mcas (items a) ), (7/V°)Bkgd (items b) ), or (7r°)Meas/(7r°)Bkgd (item c) ). The systematical errors are quoted at two pr values to give 
an indication of the dependence on transverse momentum. The errors are estimated for the narrow shower identification criterion (S2). The 
total estimated systematical error on 7Meas/7Bkgd is given as the quadratic sum of the various contributions. See Ref. [23] for full details. 
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